Hole Superconductivity in MgB2: a high Tc cuprate without Cu 
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The theory of hole superconductivity explains high temperature superconductivity in cuprates as 
driven by pairing of hole carriers in oxygen pn orbitals in the highly negatively charged Cu — O 
planes. The pairing mechanism is hole undressing and is Coulomb-interaction driven. We propose 
that the planes of B atoms in MgB2 are akin to the Cu — O planes without Cu, and that the 
recently observed high temperature superconductivity in MgB2 arises similarly from undressing of 
hole carriers in the planar boron px,y orbitals in the negatively charged B~ planes. Doping MgB2 
with electrons and with holes should mirror the behavior of underdoped and overdoped high Tc 
cuprates respectively. We discuss possible ways to achieve higher transition temperatures in boron 
compounds based on this theory. 



I. INTRODUCTION 

Superconductivity with a remarkably high transition 
temperature AQK) has recently been discovered in 
MgB2 [0. It has been proposed that it results from 
strong electron-phonon interaction and high phonon fre- 
quency due to the light ionic masses , and indeed mea- 
surement of an isotope effect has been reported |^ . Here 
we propose instead that the phenomenon is another clear 
example of the mechanism of hole superconductivity at 
work 1 1] . The theory of hole superconductivity proposes 
that this is the universal mechanism of superconductiv- 
ity in solids . 

No measurement of the Hall coefficient of MgB2 has 
yet been reported in the literature to our knowledge. The 
Hall coefficient for other metal diborides with the same 
crystal structure {YB2, TiB2, VB2, ZrB2, NbB2, etc.) 
has been measured and found to be negative . Corre- 
spondingly we expect no superconductivity in those com- 
pounds (no superconductivity in those compounds has 
ever been detected to our knowledge). Instead, we expect 
the Hall coefficient of MgB2, certainly for transport in 
the planes, to be positive, indicating the dominant role of 
hole carriers in the normal state transport expected for 
superconductors within our theory. In other words, we 
expect the family of metal diboride compounds to be an- 
other clear example of the Chapnik-Feynman empirical 
rule [^-10 1 that superconductivity exists when the nor- 



mal state transport is hole-like and does not exist when 
it is electron-like. 

In the theory under consideration here 1 , supercon- 
ductivity is an 'undressing' transition, and it can only be 
driven by carriers in bands that are almost full. It is ar- 
gued that the dressing of quasiparticles in an electronic 
energy band, due to electron-electron interactions, is an 
increasing function of the electronic band occupation. 
When the Fermi level is near the top of the band, the 
carriers (hole carriers) are most heavily dressed, and the 



normal state transport is largely incoherent; when these 
heavily dressed hole carriers pair they partially undress, 
and this is the driving force for superconductivity. The 
superconducting condensation energy is kinetic in origin, 
as undressed carriers have a lower effective mass. Sim- 
ilarly, undressing also occurs upon hole doping a nearly 
full band in the normal state. If the Fermi level is not 
close to the top of the band, carriers are not as heavily 
dressed in the normal state and do not gain enough by 
pairing to overcome Coulomb repulsion, and supercon- 
ductivity disappears. 



II. WHAT MAKES MGB2 A HIGH 
TEMPERATURE SUPERCONDUCTOR 

The crystal structure of MgB2 is the so-called AIB2 
structure: honeycomb layers of boron atoms alternate 
with hexagonal layers of Mg atoms. B atoms in different 
planes are on top of each other, and the Mg atoms are 
at the center of the hexagons defined by the B atoms. 
It was proposed long ago that in metal diborides the 
boron atoms accept electrons from the metal, and the 
boron planes become negatively charged The sim- 
plest ionic picture would suggest that in MgB2 the Alg 
donates two electrons to the B planes, and the ionic com- 
pound A/g++(_B^)2 results, which will be metallic due to 
boron band overlap. This picture is supported by band 
structure calculations. 

Within the theory of hole superconductivity what 
drives superconductivity in the cuprates is hole transport 
through the planar oxygen pir orbitals p^ . That band 
is full in the undoped cuprates, and we have proposed 
that when the materials are hole doped, hole carriers 
will go predominantly into that band (although transport 
through the O pa-Cu d^2_y2 orbitals may also occur |p^). 
For the electron-doped cuprates, we have proposed that 
holes in the planar O p-K orbitals are induced by electron 
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doping ]14[ | and that they become mobile when oxygen 
is removed through anneahng in a reducing atmosphere 

Band structure calculations for the doped cuprates are 
complicated due to the strong Coulomb repulsion on the 
Cu atoms, and so far they have not shown evidence that 
the Fermi level cuts an oxygen pi: band near its top. In 
contrast, band structure calculations for MgB2 are sim- 
pler and expected to be accurate, and several calcula- 
tions have been reported in the literature, using different 
methodologies, that are in essential agreement p6|-p^,p| . 

Figs. 1 and 2 show band structure calculation results 
of Armstrong et al for ^ and AIB2 [|l9| respec- 

tively. The essential fact that makes MgB2 a supercon- 
ductor, according to our theory, is that the Fermi level 
at the r point cuts two bands right near the top, as seen 
in Fig. 1 (bands 3 and 4 in the terminology of Arm- 
strong et al [0). These two bands give rise to nearly 
cylindrical hole-like Fermi surfaces around the F point 
(see also Fig. 3 of Ref. 2), indicating that the transport 
from these hole carriers will be dominantly in-plane. The 
atomic states giving rise to these bands are dominantly 
boron planar p^^y orbitals §,|6|. Within the theory of 
hole superconductivity, it is those heavily dressed hole 
Px.y carriers that give rise to the high temperature su- 
perconductivity of MgB2. As emphasized by Armstrong 
et al, the main difference between MgB2 and AIB2 is 
that in the latter compound these p^^y orbitals are com- 
pletely filled and do not contribute to conduction. The 
conduction in AIB2 is mainly electron-like, as suggested 
by Fig. 2 and by the measured negative Hall coefficient 
of the similar compound YB2 f^. Hence AIB2 is not a 
superconductor. 

The existence of a hole-like Fermi surface is a neces- 
sary condition for superconductivity within our model. 
However to obtain high temperature superconductivity 
the presence of negatively charged conducting structures 
is also found to be favorable, as negative ions give rise 
to larger values of the parameter At that drives super- 
conductivity 1^^. Thus the fact that the boron planes 
in MgB2 are essentially B~ , i.e. are highly negatively 
charged, naturally fits within this picture. The same sit- 
uation is found in the high Tc oxides, where the Cu — O 
planes are also highly negatively charged, with two extra 
negative charges per unit cell (Cm"'""'" (0^)2). 

Figure 3 shows the generic behavior of Tc versus hole 
concentration predicted by our theory . The magni- 
tude of the parameter At used, that drives superconduc- 
tivity in our model, is At ^ 0.37eV for the case of Fig. 3, 
well within the range obtained from first-principles cal- 
culations of this parameter for p-orbitals and negatively 
charged ions ||2^. According to the MuUiken population 
analysis of Armstrong and Perkins [|l6| the hole occu- 
pation in the boron p^^y orbitals in MgB2 is 0.07. Al- 
though the position of the maximum Tc in our model 
can vary somewhat with parameters, it generally occurs 
for lower hole concentrations, as seen in the example in 
Fig. 3. This would imply that MgB2 is slightly over- 



doped, so that a small decrease in the hole carrier con- 
centration would increase Tc up to a maximum of around 
50K . Doping with electrons, for example in a compound 
Mgi-xAlxB2 (assuming it forms with the same struc- 
ture) would bring Tc over the maximum and drive it to 
zero in the underdoped regime. Doping with holes, for 
example with Mgi^yNayB2 (assuming it forms with the 
same structure) would drive Tc to zero in the overdoped 
regime. 

The theory of hole superconductivity predicts a cross- 
over from strong to weak coupling regimes as the hole 
concentration increases a scenario which is in qual- 
itative agreement with observations in the cuprates. We 
may expect to see a similar scenario (although less pro- 
nounced) in the doped magnesium diborides. In the 
underdoped regime [M gi-xAlxB2) increasingly incoher- 
ent transport and higher resistivity, decreasing coherence 
length and increasing gap ratio as x increases, possi- 
bly even pseudogap behavior and charge inhomogeneity. 
This should only occur for a small range of x however, af- 
ter which the Px.y bands will become full and the behavior 
will change sharply for larger x: superconductivity will 
dissappear, the Hall coefficient will become negative, the 
transport will become coherent, and the resistivity will 
decrease as x increases further with increasing electron 
carriers. The lattice stability should also increase in this 
regime of large x. In contrast, in the overdoped regime 
{M gi-yN ayB2) we expect increasingly coherent behav- 
ior as y increases, with lower resistivity, increasing super- 
conducting coherence length diverging as Tc approaches 
zero, and gap ratio close to the BCS weak coupling value. 
The London penetration depth should decrease mono- 
tonically as the hole concentration increases from under- 
doped to overdoped p^ . 

III. EFFECT OF PRESSURE 

Within our theory a decreasing B — B intraplane dis- 
tance should increase Tc 0, as seen in Fig. 3. Similarly 
we have argued that in high Tc oxides the intrinsic effect 
of increase of Tc with pressure is caused by decreasing 
intraplane O — O distance . This is because supercon- 
ductivity is driven by the correlated hopping parameter 
At that depends exponentially on interatomic distance. 
Hence we expect application of pressure in the plane di- 
rection in MgB2 will increase Tc- Hydrostatic pressure 
should also increase Tc assuming it leads to a decrease 
in the B — B intraplane distance. However the situation 
could be more complicated if charge transfer between the 
B planes and the metal occurs under pressure, in which 
case the sign and magnitude of the change in Tc would 
depend on the sign of the charge transfer, whether the 
system is in the overdoped or underdoped regime, and 
on the relative weight of the change induced by charge 
transfer and change induced by changing interatomic dis- 
tances. All of these effects could be disentangled by mea- 
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suring changes in lattice constants and Hall coefficient 
under pressure. 



IV. HOW TO ACHIEVE HIGHER Tc'S IN THE 
DIBORIDES 

Various calculations of band structures in diborides 
suggest that a rigid band picture works reasonably well 
p6|- Jl9| , p4[ . In the transition metal diborides it is found 
that the Fermi level states are dominantly of metal 3d 
character [|2^, in contrast to the main group diborides 
where the boron 2px^y orbitals dominate the Fermi level 
states. According to the theory discussed here, the lat- 
ter situation is the favorable one for superconductivity, 
and the position of the Fermi level in MgB2 is close to 
optimal. Hence what remains to be optimized is the in- 
teratomic distances, which may be achieved by 'chemical 
pressure'. 

Consider the compound BeB2 |^,^. Its structure is 
similar although not identical to that of MgB2, however 
the interatomic distances should be significantly smaller 
than those in MgB2 due to the smaller size of the Be 
atom. The lattice constants for BeB2 obtained by aver- 
aging experimental data over a larger unit cell have been 
inferred to be a = 2MA, c = 2.87 A while those of 
MgB2 are a = 3.084A, c = 3.522A H]. Hence the crucial 
B — B intra-plane distance should be about 5% shorter 
in BeB2, allowing in principle for considerably higher Tc- 

However even though Be and Mg are both nominally 
divalent the ionization potentials of Be are significantly 
larger: the first and second ionization potentials are 
Ei = 9.32eV, En = 18.21ey for Be, and Ej = 7MeV, 
En = 15.03ey for Mg (28). Hence the charge transfer 
from Be to B in BeB2 will be less than that from Mg to 
B in MgB2, and the Fermi level in BeB2 will be lower 
than shown in Fig. 1, probably beyond the regime where 
superconductivity occurs. Hence the system BeB2 needs 
to be doped with electrons, and we suggest to achieve 
that through the compound Bei^xAlxB2. If this com- 
pound forms in the AIB2 structure, we predict that with 
increasing x the Tc versus rih curve of Fig. 3 will be 
mapped from overdoped to underdoped, with a maxi- 
mum Tc significantly larger than in the MgB2 structure. 

Similarly a LiB2 or NaB2 compound with the same 
structure and doped with electrons, e.g. Lii-xAlxB2, 
would map the Tc versus rih curve, for larger values of x 
[x > 1/2). The Li case should yield higher T^'s than the 
Na case due to the smaller interatomic distances. 

Finally we comment on the transition metal di- 
borides. Several compounds of the form MB2 
form with the AIB2 structure, e^g. M — 
Sc, Ti, V, Cr, Y, Zr, Nb, Mo, Hf, Ta, W [|mL AU are re- 
ported to have negative Hall coefficient [^,0 and the 
transport is dominated by carriers in the metal d orbitals 
1 18,^. The band structure calculations of Armstrong 
1 24 1 show that the Fermi level is above the Fs point in 



Figure 1 , and approaches that point as one moves to the 
right and down in the periodic table. Doping these ma- 
terials with holes, for example by forming intermetallic 
compounds where transition metals are partially substi- 
tuted by alkali or alkali earths, should bring the Fermi 
level down to the B p^.y bands and give rise to supercon- 
ductivity. However these cases should be like the two- 
band situation of ref. |^3|, where the hole carriers drive 
the system superconducting in the presence of a large 
number of electron carriers, and Tc's should be consider- 
ably lower than in AIgB2. 



V. ISOTOPE EFFECT 

The presence of an isotope effect is usually considered 
evidence that electron-phonon coupling drives supercon- 
ductivity. However, an isotope effect is also expected 
in the mechanism considered here. The parameter that 
drives superconductivity. At, depends sensitively on in- 
teratomic distances, as does the single particle hopping t 
Q . Let us assume we have 



At{q) =At + -f'q 



(la) 



(lb) 



where q is an optic phonon mode, and 7' = ^At/t. The 
zero-point motion of the atoms will yield a small increase 
in the effective At, that depends on the ionic mass: 



J 



At 



Atcff = V< Ai >2 = At H- 4- < >= At + -/^— <q^ > 



At 



t^ 



(2) 



Using < q^ >= hw/2K, with lo the phonon frequency 
and K the force constant, we find that if the ionic mass 
M changes by 5M , the parameter Atc/f changes by 



^(Ate//) 



(3) 



Finally, with D = 2zt the bandwidth {z the number of 
nearest neighbors) and A = j^jKD the dimcnsionless 
electron-phonon coupling we have 



AthujdM 

6(Atcff) = zX 

^ "^^^ t 2 M 



and the change in Tc is given by 
dlnTr dlnTr 



zXhuj 



dlnM dlnAt 2t 



(4) 



(5) 



For the parameters used in Figure 3, dlnTc/dlnAt — 
14.5, t ~ 0.63ey, and using u; = 750K as an estimate of 
the phonon frequency we find for the isotope coeffi- 
cient 
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dlnTc 
dlnM 



2.2A 



(6) 



so that to reproduce the measured isotope shift as ^ 
0.26 (I would require a rather small A, A ^ 0.12. The 
estimated value of A in ref. (2| is much larger, A ~ 0.7. 

While the calculation discussed here is not expected to 
be very accurate, it illustrates the fact that the mecha- 
nism of hole superconductivity is very sensitive to ionic 
vibrations and generically leads to a positive isotope ef- 
fect due to the increased effective At from larger zero- 
point vibrations. Eq. (6) and the measured isotope shift 
suggest that the actual A in the material may be substan- 
tially smaller than estimated in ref. Alternatively, 
the isotope shift estimate Eq. (6) may be reduced by the 
fact that a smaller ionic mass would also be expected to 
lead to a slightly larger average interatomic distance, an 
effect which would partially compensate the increase in 
effective Ai estimated here. 



VI. CONCLUSIONS 

We have proposed here an interpretation of the ob- 
served high temperature superconductivity in MgB2 
based on the theory of hole superconductivity. According 
to this theory, the essence of high cuprates is metal- 
lic oxygen, more specifically holes conducting through 
nearly full oxygen pir planar orbitals in the highly neg- 
atively charged Cu — O planes Similarly, holes 
in nearly filled boron planar px,y orbitals in the nega- 
tively charged planes should drive superconductiv- 
ity in MgB2- In this sense, MgB2 is then a beauti- 
ful realization of the essential physics of superconductiv- 
ity in cuprates, without the complications of Cu dx2_y2 
orbitals, antiferromagnetism, chains, etc., that we have 
argued obscure the physics and are non-essential in the 
cuprates [^,^ . Already in our early work we noted in con- 
nection with high Tc cuprates "a moment on the cation 
(like Cu^+) is not needed, in particular a closed shell 
should do" , and "The high Tc oxide structures with Mg 
in place of Cu, if they formed, would be excellent candi- 
dates to support our picture and rule out other mecha- 
nisms" . Remarkably, the band structure calculations 
in the borides and their interpretation P,p^ directly sup- 
port the scenario expected in our theory to be favorable 
for hole superconductivity. 

We expect that the various effects predicted by the 
theory of hole superconductivity for the cuprates [||J|] 
should exist, and be of appreciable magnitude in MgB2, 
due to its high Tc, in particular: (1) tunneling asymme- 
try of universal sign, i.e. larger current for negatively 
biased sample; (2) charge imbalance of quasiparticle ex- 
citations (quasiparticles have a net positive charge), re- 
sulting in particular in positive thermopower for NIS and 
SIS tunnel junctions; (3) apparent violation of low energy 
conductivity sum rule, i.e. a larger ^-function weight 
(smaller London penetration depth) than expected from 



the missing area in the low frequency conductivity; (4) 
color change, i.e. transfer of optical spectral weight from 
high frequencies, up to the visible range, down to low fre- 
quencies upon onset of superconductivity; (5) evidence 
for 'undressing' in angle- resolved photoemission p9| , i.e. 
emergence of a sharp quasiparticle peak in the super- 
conducting state from a weaker peak and an incoherent 
background in the normal state, as well as evidence of in- 
creased coherence upon hole doping in the normal state. 

Furthermore, as discussed in previous sections, we ex- 
pect Tc to be an increasing function of pressure applied 
in the plane direction, positive Hall coefficient for MgB2, 
and Tc versus hole concentration of the generic form 
given by Figure 3, with strong coupling to weak cou- 
pling crossover in various properties as function of in- 
creasing hole doping in the various diboride compounds 
discussed in the text. As we noted in earlier work |]l^ , ^ , 
the same generic Tc versus hole concentration behavior 
is seen in transition metal alloys, a phenomenon known 
as "Matthias' rules" ||. 

Of course there will be other theories proposed to ex- 
plain the superconductivity of MgB2 and related com- 
pounds, and in particular an electron-phonon one has 
already been advanced Even though an isotope ef- 
fect may appear to favor such a theory, we have pointed 
out here that such an effect is also expected within the 
theory of hole superconductivity. Fortunately all theories 
should not yield identical predictions as to which other 
diborides should give rise to high temperature supercon- 
ductivity, nor about what their properties are. Spelling 
out such predictions clearly, preferably before the exper- 
iments are performed, should make it easier to ascertain 
the relative merits of the various theories. 

In closing we remark that the theory of superconduc- 
tivity through hole undressing discussed here predicts the 
symmetry of the superconducting state for all supercon- 
ductors to be s-wave. We suggest that the apparent evi- 
dence for d-wave superconductivity in the cuprates may 
be connected to the presence of Cu d orbitals near the 
Fermi level. Since no such orbitals exist in MgB2 we 
expect that clear evidence for s-wave superconductivity 
will be experimentally found. 
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FIG. 2. Band structure calculation results for ^4^52 in 
the plane directions reproduced from Ref. 19 by Armstrong 
et al. In contrast to Figure 1, the hole pockets at the F 
point no longer exist as the Fs point is considerably below the 
Fermi level. The Fermi surface for this case is predominantly 
electron- like. 



FIG. 3. Tc versus hole concentration in the model 
of hole superconductivity for a two-dimensional case. This 
behavior is generic for this model. Values for the bandwidth, 
correlated hopping parameter, on-site and nearest neighbor 
repulsion used were D = 5eV, At — 0.3725el/, U = 5eV, 
V = respectively. The dashed line indicates the behavior 
expected under application of physical or chemical pressure, 
with the parameter At increased to 0.375ey. 



FIG. 1. Results of band structure calculation for MgB2 in 
the plane directions reproduced from Ref. 16 by Armstrong 
et al. Note the hole pockets at the Fs point (where the arrow 
points). They give rise to cylindrical hole Fermi surfaces in 
the c direction for the third and fourth bands. 
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